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Abstract. Experimental highlights of the RHIC program are reviewed. 

INTRODUCTION 

Relativistic heavy-ion collisions provide a unique environment to study matter under 
extreme conditions of high temperature and energy density. RHIC (Relativistic Heavy 
Ion Collider) at the Brookhaven National Laboratory in Upton, NY provides us with 
such collisions. Since the beginning of the RHIC operation in the year 2000, RHIC has 
provided us not only with Au+Au collisions at the top energy of y/s NN = 200 GeV, but 
also at 62 GeV. RHIC has also provided collisions of lighter systems such as Cu+Cu at 
y/s NN = 200, 62, and 22.5 GeV, and finally the p + p and d+Au collisions at y/s NN = 200 
GeV that were used as references. 

The amount of data obtained so far is as overwhelming as the results. I will try to 
summarize where we are, with RHIC RUNV just about to start. 

ELLIPTIC FLOW 

In non-central collisions the initial spatial anisotropy is transformed into an anisotropy 
in momentum-space if sufficient interactions occur among the constituents within the 
system. Once the system has expanded enough to quench the spatial anisotropy, further 
development of momentum anisotropy ceases. This self-quenching process happens 
quickly, so elliptic flow is primarily sensitive to the early stages of the collisions yj]. 



Hydrodynamics 

The elliptic flow V2 as a function of p t for K^,A, 0,2, and £1 is depicted in Fig. [l] 
|0,|3l]. The 0,2, and Q have low hadronic cross-sections, therefore the large V2 observed 
suggest that the elliptic flow is built up in the partonic stage. The expected range of V2 
from hydrodynamic calculations is also shown in Fig. [l] A more detailed comparison 
can be seen in Fig. |2 (left panel), where the mass dependence hydrodynamic results JU] 
are compared to the V2 measurements of n,K,p, and A ||4|]. Hydrodynamics describes 
well the mass dependence observed in the data that is characteristic of a common flow 
velocity. Since an ideal hydrodynamic fluid is a thermalized system with a zero mean 
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FIGURE 1. Azimuthal anisotropy for strange (left panel) and multi-strange (right panel) hadrons in 
minimum bias Au+Au collisions Q|. Data measured by STAR. The dashed lines show a common fit to 
the and A + A data [5]. The shaded areas are hydrodynamic calculations Jit 



free path that yields to the maximum possible vj, the good agreement between the 
measured V2 and the hydro results |Q] suggests thermalization in heavy-ion collisions 
atRHIC. 



Constituent Quark Scaling 

While hydrodynamic calculations keep increasing as a function of pj, the measured 
V2 saturate at pj > 2 GeV/c ||2j]. The saturation value for mesons is about 2/3 of that for 
barypns. This separation pattern holds for n,K,A, and E, and seems to hold for and 
£1 lylla]. This result and the baryon-meson splitting of the high pj suppression pattern 
|7|] suggest the relevance of the constituent quark degrees of freedom in the intermediate 
pj region [8]. v% scaled by the number of valence quarks n as a function of pr/n is 
depicted in Fig.|2](right panel). The lower panel of Fig.|2](right panel) displays the ratio 
between the measurements and a polynomial fit to all the data. At low pr/n the observed 
deviations from the fit follow a mass ordering which is expected from hydrodynamics. At 
higher pr, all vijn measurements are reasonably close to unity showing the constituent 
quark scaling. 

FREEZE-OUT PROPERTIES 

The measured particle spectra and yields H] and event-by-event (pr) fluctuations ifioll 
indicate a nearly chemically and kinetically equilibrated system at the final freeze-out 
stage. 



Polynomial Fit 




Transverse momentum p t (GeV/c) p/n (GeV/c) 



FIGURE 2. Left panel: V2 as a function of p t from minimum bias Au+Au collisions |4] measured by 
STAR. The dotted-dashed lines are hydrodynamic calculations using an equation of estate (EOS) with a 
first-order hadron quark-gluon plasma phase transition [ 1 ] . The description of the data worsen if a hadronic 
EOS is used 1 1]. Right panel: Measurements by STAR of the scaled V2(pr /n)/n for identified hadrons 
(upper panel) and the ratio between the measurements and a polynomial fit through all data points (lower 
panel) except the pions for ^/s^ = 200 GeV minimum bias Au+Au collisions (jj. 



Chemical and Kinetic Freeze-out Parameters 



STAR has now measured hadron distributions at */s NN = 200 and 62 GeV 111 



1211 . Chemical freeze-out properties were extracted from stable particle ratios within 
the thermal model 11511 . The extracted chemical freeze-out temperature T c \ x and the 
strangeness suppression factor y s are depicted in Fig. |3j Kinetic freeze-out properties 
from particle pj distributions were extracted within the blast wave model II 1611 . Figure 
|4| displays the extracted kinetic freeze-out temperature Tu n and the average radial flow 
velocity (j8r). The results at y/s NN = 62 GeV are found to be qualitatively the same as 
those obtained at y/s NN = 200 GeV and resonance decays are found to have no significant 
effect on the extract kinetic freeze-out parameters lfl2ll . 

T c h is independent of centrality. T^ n obtained from ft, AT, and p decreases as a function 
of centrality, while the corresponding (j8j) increases. This is evidence that the system 
expands between chemical and kinetic freeze-outs, which brings the system to a lower 
temperature. The constant T c h suggests that hadronic scatterings from hadronization to 
chemical freeze-out may be negligible because they would result in a dropping T c \ x . 

y s increases from p + p to peripheral and central Au+Au collisions. In central Au+Au 
collisions, % is ~1 suggesting that strangeness is saturated. The T c f, of 0,E, and Q. is 
higher than that of TZ,K, and p, while the (/3r) is lower. Noting that the 0,S, and Q. 
have small hadronic cross-sections, they may chemically and kinetically freeze-out at 
the same time. 
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FIGURE 3. Left panel: Extracted chemical freeze-out temperature from stable particle ratios ifTTl 
O .Right panel: Extracted strangeness suppression factor from stable particle ratios fill O . STAR 
measurements. 
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FIGURE 4. Extracted kinetic freeze-out temperature as a function of the average flow velocity within 
the blast-wave model STAR measurements. 



SYSTEM SIZE DEPENDENCE 

The comparison between the charged hadron pseudorapidity distributions measured 
by PHOBOS [13] in Cu+Cu and Au+Au with similar number of participants N par t is 
presented in Fig.|5j Both distributions are comparable within errors, showing that bulk 
particle production depends mainly on the number of participating nucleons. The same 
is true for different N part and also at ^/s^ = 62 GeV 11311 . Many other measurements 
from BRAHMS, PHENIX, PHOBOS, and STAR support this argument JUHKHS]. 
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FIGURE 5. Charged hadron pseudorapidity measured by PHOBOS for Cu+Cu and Au+Au collisions 
at y/s NN = 200 GeV with similar number of participants N pan 



LOW P T HADRONS 

The nuclear modification factor Raa of ft and p measured at y = 3.2 by BRAHMS 
and y = by PHENIX in central Au+Au collisions does not depend on rapidity lEoll . 
suggesting that the same mechanisms is responsible for the nuclear modifications. It has 
been predicted W21I1 that the magnitude of jet quenching should depend on both the size 
and the density of the created absorbing medium. The averaged pion Raa as a function 
of Npart for both forward rapidity and mid-rapidity are shown in Fig. |6j The average 
was performed in the interval 2 < pj < 3 GeV/c. The mid- and forward rapidity pion 
suppression for the most central Au+Au collisions are found to be the same magnitude. 
However, the Raa measured in forward rapidity shows significantly stronger rise towards 
peripheral collisions as compared to Raa at mi-rapidity, differing on the level of 35% for 
(N part ) ~ 100- This is consistent with the model of parton energy loss in a strongly 
absorbing medium lE2Ll23"ll . In this picture, at mid-rapidity, the emission is dominated 
by the emission from the surface which quenches the dependence of Raa on the system. 
On the other hand, at forward rapidities, the transition from surface to volume emission 
can occur, which leads to a stronger dependence on the number of participants. 



INTERMEDIATE P T HADRONS 

Intermediate pj (pj < 6 GeV/c) protons behave differently than mesons at heavy-ion 
collisions [6]. This behavior can be seen in Fig.|7](reft panel) that depicts the ratio of p to 
it spectra measured by STAR in both central Au+Au collisions and in p + p collisions. 
The large enhancement of the pjn ratio at intermediate pj in Au+Au collisions indicates 
that jet fragmentation in vacuum is not the dominant source of particle production in 
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FIGURE 6. Averaged R M in the interval 2.0 < p T < 3.0 GeV/c at mid-rapidity (PHENIX) and at 
forward rapidity (BRAHMS) as a function of number of participants at .Js NN = 200 GeV fj2j 



this pr range. An enhancement, roughly peaked at at the same position, is also observed 
in Rep, as shown in Fig. |7] (right panel). Here, Rep is the ratio between central and 
peripheral Au+Au collisions normalized to the binary collision scaling expectation. 
Like V2 [2], Rcp also separates into baryons and mesons. The K* LI 7], which is a 
meson with the mass close to the proton mass, follow the behavior of mesons. This 
proves that this separation is not due to the hadron mass. This grouping is violated in 
RaAi where the reference is from p + p collisions rather than peripheral collisions. A 
strong enhancement in strange baryon Raa with increasing enhancement for increasing 
strangeness content is observed at intermediate pr H 1 811 - This enhancement is depicted 
in Fig. [HI 



HIGH P T HADRONS 

The suppression of high pj hadrons in central Au + Au collisions was one of the 
unexpected and important phenomena observed at RHIC. The Raa of 7T° measured by 
PHENIX in central Au+Au collisions at ^/s^ = 200 GeV Q is presented in Fig. 
The suppression is quite strong (Raa ~0.2) and remains approximately flat up to 20 
GeV/c. Partonic radiative energy loss models lEM I26L l27ll reproduce this behavior well. 
In particular, one calculation depicted in Fig. |9] uses an average gluon average density 
dPL/dy ~ 1200 B. The Raa measured by PHENIX and STAR are shown in Fig. 
II24II and Fig. [10] LI 9], respectively. The difference between charged hadrons Raa an d the 
K Raa in the intermediate pr region (pj < 6 GeV/c) is due to the proton contribution 
as observed in Fig. [I0|and nicely described by recombination models. This behavior 
disappears for p T > 6 GeV/c and are also explained by the same recombination models 
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FIGURE 7. Left panel: p/n ratio measured by STAR as a function of pj for central Au+Au (0-5%), 
af+Au and p + p collisions 1 19]. For the Au+Au measurements, the error bars are statistical only and the 
solid lines are the systematic uncertainties. For the d+Au and p + p measurements, the errors shown are 
combined statistical and systematic. Right panel: Rep as a function of pj for identified particles measured 
by STAR fl8lfT^l . Errors are statistical and systematic. Grey bands are common scale uncertainties from 
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statistical and systematic. Grey bands are common scale uncertainties from Nonary 



JET CORRELATIONS 



A well defined back-to-back peak that is characteristic of di-jets is observed with negligi- 
ble background in Fig.[TT]in both peripheral and central Au+Au collisions measured by 
STAR IB ill . While the yield is substantially less in central Au+Au collisions, the widths 
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FIGURE 9. Raa as a function of pr of 7T° and charged hadrons measured by PHENIX in central Au+Au 
collisions (0-10%) with theoretical predictions 1 25, 26]. The shaded areas are the systematic uncertainties. 
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FIGURE 10. Rcp for the 5% most central collisions measured by STAR, normalized by the peripheral 
60-80% collisions, as a function of pj 1 1 911 - The bands show combined statistical and systematic errors. 



of the back-to-back peaks appear to be independent of centrality. The away-side hadron 
triggered fragmentation functions [26] for d+Au and Au+Au collisions IBlll measured 
by STAR as a function of ZT = pf soc / p f g is depicted in Fig.[j2](left panel). Fig.EKleft 
panel) shows that the shape of the away-side fragmentation functions is unchanged from 
J+Au to central Au+Au collisions. However, the yields are reduced by a factor of ~4 
in central Au+Au collisions. Even though the shape is consistent with prev ious predic- 
tions in the zt range measured, the magnitude is smaller than expected Il2fl1 . A different 
calculation predicted that significant energy loss should be associated with significant 
broadening of the away-side hadron azimuthal distribution [13211 . in contradiction to the 
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FIGURE 11. Azimuthal distributions of away-side charged hadrons for 8 < p 1 j S < 15 GeV/c in 
two different centralities at y/s NN = 200 GeV Au+Au collisions. These are the total distributions, no 
backgrounds have been subtracted. 



STAR measurements. 

When low-pj associated particles are observed in central Au+Au collisions opposite a 
high-pr trigger particle, they are found to be significantly broadened in and enhanced 
in number compared to p + p collisions 113311 . Contrary to the maximum found in p + p 
and d+Au collisions, STAR measured the (pj) of the associated particles appears to 
be a minimum at A0 = % in central Au+Au collisions l34ll . In addition, the away- 
side associated particle yield is flat or may have a small dip at A0 = n Q4I1 . These 
phenomena are displayed in Fig. [12] (right panel) and have led to predictions that we 
may be observing jets that have been deflected by radial flow or a Mach cone effect 
associated with conical shock waves ll35ll . 



HEAVY FLAVOR 

The nuclear modification factor Raa of non-photonic electrons measured by PHENIX 
ll36ll and STAR 0] in central Au+Au collisions is shown in Fig. [13] Both experiments 
observe a very strong suppression of the non-photonic electrons. The suppression has 
approximately the same shape and magnitude as the suppression for hadrons, which is a 
quite surprising result, since the massive quarks are expected to radiate much less energy 
than the lighter u and d quarks. It is important to confirm and complement these results 
with direct measurements of open charm and open bottom contributions. The first direct 
reconstruction of D mesons by STAR IB7ll is the first step in this direction, and more 
improvement will come with the implementation of new detectors from both PHENIX 
and STAR. 

The V2 of non-photonic electrons as a function of px ll3^1 measured by PHENIX is 
shown in Fig. [14] (left panel). The non-zero flow measured at pj < 2 GeV/c, where the 
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FIGURE 12. Left panel: Top: Hadron triggered fragmentation functions dN/dzr for away-side charged 



hadrons as a function of zt for 8 < p T < 15 GeV/c at ^/s^ = 200 GeV d+Au and Au+Au collisions. 
The solid line is an exponential fit to the zt distribution for d+Au. The dashed lines represent the same 
exponential fit scaled down by factors of 0.54 and 0.25 to approximate the yields in 20-40% and 0-5% 
Au+Au collisions. Bottom: Ratio of the hadron triggered fragmentation functions for Au+AuAf +Au. Right 
panel: Azimuthal distributions of associated charged hadrons for two different centralities at y/s^ = 200 



GeV Au+Au collisions. The histograms indicate the systematic uncertainty bands. 







f=S systematic error 




Au+Au @),>fs = 200 GeV, 0-10% Centrality 


statistical &n 




< 


- PHENIX PRELIMINARY 


I n - " " 1 11. 


1 Tj A 






uncertainty 


i p*pref. 


1.2 


~ (1)q_hat = 0GeV 2 /frn 





1 

0.B 



0.2 







(4) dN a / dy = 1000 












(2) q hat » 4 GeV 2 /fm 






(3) q_hat= 14 GeV 2 /fm 


"l I I I I 


H 1 1 I 1 1 1 I h r i i I ■ i i i 


I. I,. I, I MM 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

P T IGeV/c] 



< 3r 
< 

rr 

2.5 



1.5 



0.5 



• STAR d+Au 

» STAR Au+Au (0-5%) 

STAR preliminary 



iil 



2 



8 9 10 

P T (GeV/c) 



FIGURE 13. Left panel: R/^ of non-photonic electrons measured by PHENIX 1 36] for central Au+Au 
collisions at y/s NN = 200 GeV. Right panel: Raa of non-photonic electrons measured by STAR |7] for 
central Au+Au collisions at 



^ = 200 GeV. 




FIGURE 14. Left panel: vi of non-photonic electrons attributted to semi-leptonic open charm decays 
measured by PHENIX 1 36] and compared to theoretical predictions 1 38]. Right panel: Compilation of the 
total charm cross-section production per nucleon collision a^ d as a function of energy compared to the 
NLO pQCD calculations iH. 



yield is dominated by semi-leptonic decays of open charm, suggests a sizeable flow of 
D mesons. The comparison to calculations with and without charm flow [38] shown in 
Fig .PHI (left panel) favors the interpretation of charm flow. The consequence is a strong 
interaction with the medium and a high degree of thermalization of c quarks, favoring 
the strongly coupled QGP. 

PHENIX has calculated the non-photonic electron measurements in p + p and min- 
imum Au+Au collisions at ^JJ^ = 200 GeV into a total cross-section of charm pro- 
duction per nucleon collision o^ N 113911 . STAR has calculated the combined electron 
and direct D meson measurements in minimum bias d+Au and minimum bias Au+Au 
collisions [37]. There is a factor of ~2 difference in their minimum bias Au+Au cross- 
sections. It is important to notice that STAR is sensitive to 80% of the total charm cross- 
section, while PHENIX is only sensitive to 15%. A compilation of charm cross-sections 
o^ N as a function of collision energy is displayed in Fig. [14] (right panel). Results from 
PHENIX and STAR are also shown. 



J/y 

PHENIX measures J/\y in p + p, <i+Au, Cu+Cu, and Au+Au collisions at y/s NN = 200 
GeV and Cu+Cu at y/s NN = 62 GeV at mid-pseudorapidity \ f] \ < 0.35 through the e + e~ 
decay channel and forward pseudorapidity \t]\ E [1.2,2.2] through the decay 
channel II4011 . The Raa as a function of centrality N part for the measurements at y/s m 
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FIGURE 15. Raa of J / y/ as a function of N part in d+Au, Cu+Cu, and Au+Au collisions at y/s NN = 
200 GeV measured by PHENIX. L eft p anel: The measurements are compared to models that explain the 
J/y/ NA50 anomalous suppression 14 2ll43l 14411. Right panel: The measurements are compared to models 
involving either J/y/ regeneration by quark recombination 14 3i l44l l45l 14611 or J/y/ transport in medium 
0. 



= 200 GeV are depicted in Fig.[T5j The suppression is clear and within errors it is seems 
to be independent of the collision system. Furthermore, the magnitude of the suppression 
of ~3 for the most central collisions is similar to the suppression observed at SPS ll4lll . 

The left panel of Fig. [15] shows that models that were able to explain the anomalous 
J/y/ suppression at the SPS and that were based on interactions with comovers [42], 
color screening 14311 . or QCD-inspired in-medium effects lf44fl predict a stronger sup- 
pression at RHIC. 

The discrepancy seems to be resolved by invoking the regeneration of J/y at the 
later stage of the collision via recombination of c and c quarks, which are produced 
more abundantly at RHIC. Several attempts that combine suppression and recombination 
ll43l E4I |45[ l46h reproduce the data reasonably well, as shown in the right panel of Fig. 

El 



LOW-MASS DI-ELECTRON 

The left panel from Fig. [16] depicts the measured di-electrons pairs, the background, 
and the subtracted spectra with uncertainty from PHENIX [48]. The right panel of Fig. 
[16] shows the data compared to the hadronic cocktail 148J, and to theoretical predic- 
tions [49], where the invariant mass spectra have been calculated using different 
in-medium p spectral functions and an expanding thermal fireball model. It is still pre- 
mature for any statement, but it is definitely promising. 
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FIGURE 16. Left panel: The unlike sign mass spectrum together with the subtracted spectrum. Right 
panel: The data compared to the hadronic cocktail |48], and to theoretical predictions 1 49], where a p 
spectral function is introduced with and without in-medium modifications. 



SUMMARY 

The baryon-meson effect observed at RHIC favors the constituent quark coalescence 
as the production mechanism at intermediate pj, suggesting a partonic state prior to 
hadronization. 0,2, and Q. have small hadronic cross-sections, therefore the large V2 
observed suggest that it is built up in the partonic stage. Furthermore, the thermal 
parameters extracted from these particles suggest that they chemically and kinetically 
freeze-out at the same time. 

Intermediate pj (pr < 6 GeV/c) protons behave differently than mesons at heavy- 
ion collisions. The large enhancement of the pjn ratio at intermediate pj in Au+Au 
collisions indicates that jet fragmentation in vacuum is not the dominant source of 
particle production in this pr range. 

The suppression of high pr hadrons in central Au +Au collisions was one of the 
unexpected and important phenomena observed at RHIC. The suppression is quite strong 
and remains approximately flat up to 20 GeV/c. Partonic radiative energy loss models 
11251 l26l 12711 reproduce this behavior well. 

The quality of the RHIC data is just marvellous, and the knowledge that this data has 
provided on the dense matter formed in relativistic heavy-ion collisions is without any 
doubt unexpected. We can only wonder how new data with new detector upgrades will 
surprise us. 
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